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Abstract: We study the thermodynamics of two-color QCD at high temperature and/or 
density using a dimensionally reduced superrenormalizable effective theory, formulated in 
terms of a coarse grained Wilson line. In the absence of quarks, the theory is required 
to respect the Z(2) center symmetry, while the effects of quarks of arbitrary masses and 
chemical potentials are introduced via soft Z(2) breaking operators. Perturbative matching 
of the effective theory parameters to the full theory is carried out explicitly, and it is argued 
how the new theory can be used to explore the phase diagram of two-color QCD. 
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1 Introduction 

Dimensionally reduced effective theories are known to provide a physically transparent and 
computationally efficient way to access the long distance dynamics of strongly interacting 
matter at high temperatures [1]. Their construction requires the existence of a clear scale 
hierarchy between the hard scale irT and the soft scale gT (as well as the magnetic scale 
g 2 T), where T is the temperature and g the gauge coupling (cf. the discussion and references 
in ref. [2]). At temperatures far above the deconfinement transition, this is indeed guar- 
anteed by asymptotic freedom, but there are indications that at least a modest hierarchy 
exists all the way down to the phase transition region [3] . A lot of effort has subsequently 
been devoted to the study of Electrostatic QCD (EQCD) [4-6], both perturbatively and 
using lattice simulations (see e.g. refs. [7, 8]). 

Despite the above, it is commonly believed that dimensionally reduced theories have 
little or nothing to say about the deconfinement phase transition of QCD. This is to a 
large extent due to the fact that the construction of EQCD explicitly breaks the center 
symmetry, which is known to play an important role in the phase transition dynamics even 
with dynamical fermions present [9] . For the case of pure Yang-Mills theory, this issue was 
addressed in refs. [2, 10, 11], where a new type of center-symmetric dimensionally reduced 
effective theory (ZQCD) was proposed. Such an effective theory is expected to provide a 
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quantitatively reliable description of QCD thermodynamics down to temperatures close to 
the deconfinement transition. A key ingredient in this work was the replacement of the 
temporal gauge field of EQCD by a coarse grained Wilson line, at the same time retaining 
the superrenormalizability of the theory (for closely related work, see also refs. [9, 12, 13]). 
In particular, ref. [11] demonstrated that upon perturbatively matching two-color ZQCD 
to the full theory at high temperatures, one obtains a consistent effective theory for SU(2) 
Yang- Mills theory, exhibiting — as expected — a second order confining phase transition. It 
was further found that the predictions of the theory are to a very good accuracy independent 
of the value of the single hard (order g°T) parameter that was introduced in the coarse 
graining procedure. 

The success of ZQCD as an effective theory for high temperature Yang-Mills theory 
naturally raises the question, whether similar methods can be used to build an effective 
theory for full QCD with dynamical quarks, despite the fact that the latter explicitly 
break the center symmetry. While this type of an approach would naturally not capture 
the dynamics relevant for the breaking of chiral symmetry, it should serve as a useful 
starting point for studying the thermodynamics of QCD at least when the quark masses 
are sufficiently large. Constructing and studying such a theory for the case of two-color 
QCD is the main objective of this paper. 

Alongside with the inclusion of fermions, we will in this article discuss another impor- 
tant issue related to the original construction of ZQCD in refs. [2, 11]. As we will explicitly 
show in the following sections, these original works contain an identical algebraic error that 
affects the perturbative matching of the effective theory parameters. In particular, the do- 
main wall tension and width that were previously used to determine two of the effective 
theory parameters will be shown to be uniquely determined by the perturbative matching 
of the ZQCD Lagrangian to EQCD, and thus become genuine predictions of the theory. 1 
This discovery has little significance for the qualitative findings of refs. [2, 11], yet it has 
important implications for how the matching process should be carried out. 

The paper is organized as follows. In section 2, we remind the reader of the construction 
of center-symmetric effective theories a la refs. [2, 11], after which we discuss the new Z(2) 
breaking operators that emerge upon inclusion of fundamental fermions in the case of two 
colors. In section 3, we then identify the degrees of freedom of the effective theory and 
explain how they are related to those of EQCD, while in section 4, the detailed matching of 
the ZQCD parameters to the full theory is carried out. Section 5 finally contains the first 
predictions of the new theory, followed by a concluding section 6, in which we discuss the 
prospects of using the effective theory as a tool for studying the thermodynamics of the full 
one. Some technical details are left to the appendices; as a particularly useful byproduct of 
our calculations, we will in appendix A present analytic expressions for all N c = 2 EQCD 
parameters for general values of the temperature, quark masses and chemical potentials. 

1 Whilc we show this explicitly only for gauge group SU(2), the SU(3) case follows in complete analogy. 
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2 The effective theory 



Our purpose in this paper is to formulate a three-dimensional superrenormalizable effective 
theory for two-color QCD, valid for high temperatures and finite quark chemical potentials. 
We require that in the absence of fermions the theory respect the Z(2) center symmetry, 
and that the effects of fermions can be incorporated through Z(2) breaking operators that 
are soft, i.e. suppressed by some power of the gauge coupling. To this end, we follow the 
path laid in refs. [2, 11] and introduce as our degrees of freedom a coarse grained Wilson 
line Z{x) and a magnetic gluon field A(x). The gauge transformation acts on these fields 
as 

Z(x) -> s{x)Z{x)s{x)\ A(x) s{x)[A{x) + iV]s{x)\ (2.1) 
where s{x) £ SU(2), while the Z(2) transformation acts on the Z field as 

Z(x) -> e inn Z(x) = ±Z(x). (2.2) 

The relation of the latter transformation to the SU(2) gauge invariance of the full theory 
is explained in detail in appendix B. 

As discussed at length in ref. [11], a unique feature of the gauge group SU(2) is that 
the coarse graining procedure almost preserves the group property of the Wilson line, as 
an arbitrary sum of SU(2) matrices is itself an SU(2) matrix up to a multiplicative real 
factor. This implies that we may parameterize the field Z in the form 

Z = -(eI + OT-ct), (2.3) 

where H ■ a = n a cr a (a a denoting the Pauli matrices), and £ and n a , a = 1,2,3, are 
real scalar fields. Out of these four degrees of freedom we expect three to correspond 
to the light adjoint Higgs field of EQCD, while one should be an unphysical auxiliary 
field that effectively decouples from the dynamics of the light fields and has a mass of 
the order of the cutoff scale of the effective theory (~ T), corresponding to the inverse 
length scale introduced by the coarse graining. While the heavy auxiliary field decouples 
from the dynamics in the infrared, its fluctuations in the ultraviolet render the theory 
superrenormalizable, providing important technical simplifications. 

To obtain the Lagrangian density of the effective theory, we collect all superrenor- 
malizable operators up to fourth order in the fields that respect three-dimensional gauge 
invariance. 2 This leads to the expression 
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X - Tri^g + Tr(A^A^) + V(Z) 



(2.4) 



where 53 is the effective theory gauge coupling, Di = di — i[Ai, ■ ], Fij = diAj — djAi — 
i[Ai,Aj], and the potential V(Z) reads 

v{z) = + 6 2 n 2 + ciE 4 + c 2 (n 2 ) 2 + c 3 s 2 n 2 + diE 3 + d 2 £ii 2 . (2.5) 



2 While in a three-dimensional theory operators of order five in the fields are in principle still relevant 
and those of order six marginal, we exclude them from our consideration, as their contributions to physical 
quantities are suppressed by (large) positive powers of the ratio of the effective and full theory energy scales. 
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Here, all terms with the exception of the last two operators in the potential respect the 
Z(2) center symmetry, and were present in the model constructed for pure SU(2) Yang- 
Mills theory in ref. [11]; the Z(2) violating operators then clearly result from the presence 
of quarks. It is a straightforward exercise to verify that eq. (2.4) really is the most general 
Lagrangian compatible with the required symmetries: simple redefinitions of the fields 
allow us to combine independent kinetic terms for E and II and remove a Z(2) breaking 
term linear in E, while a term cubic in II is forbidden by the vanishing of the symmetric 
structure constant d a bc in SU(2). 

Next, we follow the rationale of refs. [2, 11] and split the effective theory potential into 
'hard' and 'soft' parts, parameterized respectively by the O(g ) constants hi and Sj. This 
results in the alternative expression 

V(Z) = hi Tr(Z^Z) + h 2 (TrZ^Z) 2 + g 2 [|-fi 2 + |(fi 2 ) 2 + s 3 X 4 + yS 3 + |EI1 2 

(2.6) 

where we have assumed the Z(2) breaking couplings di to be of the soft type. The kinetic 
terms as well as the hard part of the potential possess an extended (global) SU(2)xSU(2) 
invar iance, 

z^^zn 2 , g su(2), (2.7) 

which will later be seen to translate to a shift invariance of the light physical fields of ZQCD 
upon integrating the heavy one out. As a consequence, the hard part of the potential is 
minimized by all matrices that are special unitary up to a common real factor. It is only 
the soft terms that provide the 0(g 2 ) structure inside this 'valley', necessary to match the 
effective theory potential to that of the full theory Wilson line, cf. eq. (A.l). For later use, 
we note that comparing the two sets of coupling constants produces the relations 

h = g^i) b 2 = ^(hi+glsi), 

ci = jh 2 + gj.s 3 , c 2 = j(h 2 + gls 2 ), c 3 = ^h 2 , (2.8) 
d i = ^fi s 4, d 2 = ^gls 5 . 

3 Identification of the fields 



It is well known that at high temperatures, where the renormalized gauge coupling becomes 
small, the Wilson line effectively freezes to the global minimum of its perturbative effective 
potential [14-16] (for more recent work regarding the effective potential of high-temperature 
QCD, see e.g. ref. [17]), and to correctly describe its long distance dynamics it is sufficient 
to consider only its small fluctuations around it. It is thus natural to require that the 
predictions of our effective theory reduce to those of EQCD in the same limit, as the 
Lagrangian of EQCD can be obtained from an expansion of the Wilson line potential in 
powers (and derivatives) of the temporal component of the gauge field, Aq. This property 
can be most straightforwardly ensured by explicitly integrating out the heavy degree of 
freedom in the vicinity of one of the minima of the ZQCD effective potential, and by 
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matching the resulting non-center-symmetric (even in the absence of fermions) theory to 
EQCD. Through this procedure, the light field of ZQCD becomes associated with the 
adjoint Higgs Aq of EQCD, and we automatically obtain the values of several of the effective 
theory parameters. 

In this and the following section, we will explicitly perform the high temperature 
matching of ZQCD to EQCD, and find the values of the Sj, i.e. the soft parameters of 
ZQCD. We begin this by parameterizing the field Z as in eq. (26) of ref. [11], 

Z = ^l + |(01 + if-a), (3.1) 

which amounts to the redefinition £ = v + g^(p and II = g$x, where v is a real positive num- 
ber, chosen so that (2) = (v/2)l. Clearly, the precise choice of the parametrization of Z 
can have no effect on the physics, as long as it contains the correct degrees of freedom: once 
the effective theory is matched to the full theory properly, it will automatically reproduce 
the correct long distance physics. One should nevertheless note that, had we chosen to use 
a non-linear field parametrization, we would have had to consider the Jacobian associated 
with the change of variables in the defining path integral of the theory. 

Upon rewriting eq. (3.1) as Z = ^f 3 [1 + i{gs/v)x- a] + 0{g 2 ) and comparing with 
the full theory Wilson line, 



Q(x) = V exp 



i>9 



fP 1 f/3 

/ <1tA (t,x) =t+ig dTA (r,x) + O(g 2 ), (3.2) 
J J J 

we identify the real scalar field (f> as the auxiliary heavy degree of freedom of the effective 
theory (to leading order). Subsequently, we associate the field Xa with the light, physical 
field that corresponds to the adjoint scalar Aq of EQCD, which, together with the identi- 
fication of the effective theory gauge coupling, g 2 = g 2 T + 0(g 4 ), fixes the leading order 
value of the parameter v, 

v = 2T + 0(g 2 ). (3.3) 

Beyond the field identifications, the matching of ZQCD to the full theory is performed 
by demanding that the long distance behavior of static gluonic correlators is correctly 
reproduced by the effective theory, order by order in a weak coupling expansion. As the 
effects of the soft couplings Si are suppressed by a factor g 2 in comparison with the hard 
ones hi, one- loop graphs with only hard vertices enter the effective theory calculation with 
the same power of g% as tree graphs containing one soft vertex. This implies that to obtain 
the correlators in a consistent manner, we need to determine the one-loop effective potential 
of ZQCD. This function can be read off from eq. (18) of ref. [11] for the case of pure SU(2) 
Yang- Mills theory, and it has been generalized to include the effects of fermionic operators 
in appendix C of this paper. Inspecting the result reveals the anticipated effect that the 
degeneracy of the two minima present in the center symmetric case is broken by the nonzero 
values of s^. Without lack of generality, we may choose S4 < so that the ground state 
expectation value of X is positive. Solving for the minimum of the potential iteratively, we 
find then v = vq + g 2 V2 + • • • , where 



V0 = X I~J^ , V2 = ~ i ( 4 «3«D + £*4) +-L^h~2- (3.4) 
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Comparing this to the identification made in eq. (3.3), we infer from here the first nontrivial 
relation amongst our effective theory parameters, 

h + 4T 2 /i 2 = 0. (3.5) 
4 Matching of the soft parameters 

In this section, we will perturbatively determine the values of the soft ZQCD parameters 
Si. We begin this in section 4.1 by integrating out the heavy field <fi from the effective 
theory, requiring that the resulting Lagrangian for x agrees with that of EQCD. After this, 
we will in section 4.2 match the remaining soft effective theory parameters by demanding 
that the global structure of the one-loop effective potential of ZQCD agrees with that of 
the full theory. 

4.1 Perturbative matching of the Lagrangians 

Our first goal will be to explicitly integrate out the heavy auxiliary field (j> in order to 
obtain an effective potential for x om yj to be compared with the effective potential of the 
Aq field in EQCD. At the level of the quantum (Wilsonian) effective action, integrating 
out a given field amounts to eliminating it using its equation of motion or, equivalently, 
adding to the action of the other fields all tree-level Feynman graphs containing this field 
in the internal lines and all other fields as external legs. 

At this point, the SU(2)xSU(2) invariance of the hard part of the Lagrangian (2.4) 
proves its utility. The fields x can namely be identified as the Nambu-Goldstone bosons 
stemming from the spontaneous breaking of this extended symmetry by the nonzero ex- 
pectation value (Z). As a consequence, any contribution to the static correlators of x must 
come with at least one factor of Sj, and in particular, the one-loop part of the effective 
potential of ZQCD — the second line of eq. (C.3) — need not be taken into account. In 
addition, the gauge-fixing dependent part of the effective potential matches automatically 
to EQCD. 

Keeping only terms up to fourth order in x an d rescaling the spatial gluon field to 
achieve canonical normalization of its kinetic term, we arrive at a Lagrangian for this field 
that has the exact same form as that of EQCD, 

Aight = \ Tr ^ + \{B iX f + \m\t + ^(X 2 ) 2 + • • • , (4.1) 
with the mass parameter and quartic coupling reading 

m x = 9 3 ( s i ~ 4s 3^o ~ 7} S 4 V + s 5 v 
A = 2 5 ff S2 + 4 S3 + ^i-^y 

Note that when expressed in terms of Feynman diagrams, the quartic coupling A consists 
of two contributions, one from a soft operator of the (x 2 ) 2 type, and another from a soft, 
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SU(2)xSU(2) breaking mass correction to a <f> propagator connecting two hard cubic 0% 2 
vertices. The latter contribution was missed in refs. [2, 11], which resulted in an incorrect 
matching condition for the quartic coupling. The above expression for A thus corrects the 
corresponding expression in eq. (32) of ref. [11] and further generalizes it by including the 
effects of explicit Z(2) breaking due to dynamical quarks. 3 

Finally, the expressions in eq. (4.2) can be equated with their EQCD counterparts, 
given in eq. (A. 4). This gives us two new matching conditions, 



where the constants k« , parametrizing the effects of the quarks, are defined in eq. (A. 6). 
4.2 The Z(2) breaking parameters 

The matching conditions of eq. (4.3) should be viewed as fixing the values of two linear 
combinations of si 2,3 — an interpretation that becomes trivial in the limit of unbroken 
center symmetry. In contrast, a third, independent linear combination of these parameters 
does not affect the physics of the soft scale at the leading order at all, and can thus take 
any value. This is because in the nonlinear version of our theory, where the heavy mode 
has been integrated out in a center-symmetric fashion, the three operators multiplying 
the coefficients si,2,3 are not independent, but there is a linear relation between £ , LI 2 
and (LI 2 ) 2 . As the linear and nonlinear models describe the same long distance physics, 
there must be one combination of £1,2,3 that is left undetermined by the leading order 
matching of the linear theory, and can only be found through a higher order computation. 
The insensitivity of the long distance physics to this linear combination will be further 
demonstrated in section 5, where the domain wall solution of the field equations of motion 
is discussed. 

As EQCD violates the center symmetry explicitly, it is clear that the parameters 54,5, 
which facilitate the soft breaking of this symmetry in the presence of fermions, cannot be 
found by matching to EQCD. To determine their values, we instead have to consider the 
global structure of the ZQCD effective potential, which we do by applying the Nielsen theo- 
rem [18] and concentrating on the second stationary point (local minimum) of the effective 
potential that provides additional gauge invariant observables. A natural measure of the 
center symmetry breaking is the energy density difference of the absolute and metastable 
minima, which on the full theory side is represented by the parameter 5, defined in eq. (A. 7). 
In the effective theory, the stable and metastable minima are to the leading order located 
at £ = ±fo (and n o = 0). A comparison of the values of the potentials gives then the 

3 It should be additionally noted that the same, previously missed diagram with two cubic vertices leads 
to the generation of kinetic terms of the type X 2 (DiX) 2 and (x ■ DiX) 2 with couplings of order O(gi), which 
enter the EQCD Lagrangian only at order O(gi). Such terms can in principle be canceled by adding similar 
(non-renormalizable) operators to the effective theory Lagrangian of eq. (2.4). 



Si - 4s 3 f - -S 4 V + S 5 V 



2T Tkq 




3 IT 2 ' 

2 K>2 



(4.3) 



3vr 2 T 12tt 2 T 
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matching condition s^VqT = — <5T 4 which, using eq. (3.3), leads to the identification 

s 4 = -^ = 2^2(kI 2 -k1 2 ). (4.4) 

The last parameter to be fixed, S5, does not contribute to the energy difference of the 
two vacua, but does affect the shape of the effective potential. One simple and gauge invari- 
ant (although by no means unique) quantity sensitive to S5 is the difference of the squared 
mass parameters at the two minima, i.e. m£ of eq. (4.2) and the analogous parameter at the 
metastable minimum. The corresponding quantities are also straightforward to evaluate in 
the full theory, see appendix A, yielding the last perturbative matching condition, 

2s 5 - 3s 4 = 7^2(^0 - K o )• (4-5) 

It is interesting to note that the right hand side of this equation is proportional to the 
second derivative of the parameter 5 with respect to the chemical potential(s), i.e. the 
difference of the quark number susceptibilities in the stable and metastable vacua of the 
theory. 



5 Extended field configurations 

Having now finished the leading order matching of ZQCD to the full theory, it is important 
to test its predictions in particular for quantities that are sensitive to the center symmetry. 
Perhaps the most straightforward such test is to study extended gauge field configurations, 
which probe the global structure of the effective potential. In the absence of fermions, 
and thus Z(2) breaking operators, we can construct a stable domain wall joining the two 
physically equivalent minima of the theory. With fermions, this is no longer possible, as 
the minima are not degenerate in energy, but one can still look for a rotationally invariant 
three-dimensional solution that represents a bubble of the stable vacuum in a metastable 
environment. Although this bubble evolves with time, its growth rate can be estimated 
using a semiclassical static solution, representing a stationary point of a three-dimensional 
effective action with suitable boundary conditions [19]. 

Consider first a bubble wall configuration in the full theory, in which the (static) 
temporal gauge field Aq depends only on the radial coordinate r and points in the same 
direction in color space everywhere. It is described by a single scalar function a(r), whose 
action has the form 

/•oo 

5 cfr = /3 / dr4vrr 2 
Jo 

where the potential V e ff is obtainable from eq. (A. 3). One may then solve the equation of 
motion stemming from this action with the boundary condition a(oo) = 2irT/g, and use it 
to obtain the domain wall energy density and tension as well as the bubble profile, as has 
indeed be done in refs. [20, 21] (see also ref. [22] for a recent similar calculation). 

Within the effective theory, we first note that in order to minimize the energy cost of 
creating a bubble, it is clearly optimal to have the fields S, II minimize the hard part of 
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Figure 1. Free energy profiles of the leading order domain wall solution in the center symmetric 
limit as a function of the dimensionlcss length variable x = gTr, with the wall residing at x = 0. The 
solid red curve corresponds to the prediction of the effective theory, while the dashed black one is 
the full Yang-Mills result. The boundary condition for this one-dimensional solution is a(— oo) = 
and a(+oo) = 1. 

the potential (C.3) everywhere in space, i.e. have them satisfy X 2 + I1 2 = Vq. Recalling the 
identification of eq. (3.3), we see that we can express the fields in terms of one dimensionless 
function a, ranging from to 1, as E = vq cos(7ra), | XT | = t>osin(7ra). Plugging these 
formulas into eq. (C.3), we obtain upon a trivial shift the potential 

2 4 3 

Vcs(a) =^(si - As z vl) sin 2 (7ra) + ^(s 2 + 4s 3 ) sin 4 ( TO ) - ^-| sin(^a)| 3 + 

2 3 3 C") 

+ ~^ s 4 [cos 3 (7ra) — l] + ~7j- s 5 cos(-7ra) sin 2 (7ra), 

using which the bubble profile can again be solved. 4 

Specializing for the moment to the domain wall calculation in the Z(2) invariant, pure 
Yang- Mills case, we observe that the potential of eq. (5.2) depends exactly on the two linear 
combinations of 5^2,3 that were determined in our perturbative matching, cf. eq. (4.3). We 
conclude (in contrast to the claims of refs. [2, 11]) that the domain wall tension and profile 
become genuine predictions of the effective theory. Indeed, from eq. (4.3) we infer the 
results S! - 4s 3 ug = 2T/3 and s 2 + 4s 3 = 1/(3tt 2 T), using which we straightforwardly 
obtain as the domain wall tension 

rp3 

crw 4.899 x — «0.91ctym, (5-3) 



inserting the values of the soft parameters to the domain wall potential in the N{ — case, one sees 
that the different terms combine to V c g(a) — |T 3 sin 2 (tyo) [l — i | sin(7ra) |] 2 , closely reminiscent of the full 
theory domain wall potential (see in particular the discussion in ref. [23]). We would like to thank Chris 
Korthals Altes for pointing out this issue to us. 
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where ctym = (§) 1L ~g~ denotes the full theory result [20]. In figure 1, we plot the full and 
effective theory domain wall profiles, which we find to agree at a satisfactory level. These 
curves should be contrasted with those in figure 1 of ref. [11], where it was erroneously 
argued that one can fit the domain wall tension and width to the full theory values, based 
on the algebraic error discussed in section 4.1. 

Let us finally return to the three-dimensional bubble solution, relevant when dynamical 
quarks are present. Its formation and profile are determined by a balance between a volume 
energy gain, scaling like 5i? 3 (with the energy density difference 5 introduced in eq. (A. 7)), 
and a surface energy cost, scaling like aR 2 , where R stands for the bubble radius in units 
of l/(gT). We will not attempt a full numerical solution of the corresponding equation 
of motion, which is straightforward but not particularly illuminating, but instead provide 
an analytic approximation valid in the limit of parametrically small 5, corresponding to 
weak Z(2) breaking effects. As we expect R to scale like 1/5, the explicit r-dependence 
of the action becomes then negligible and the bubble profile calculation reduces to the 
type of domain wall problem encountered above. This is usually called the thin wall 
approximation [24]. 

In the thin wall approximation, the bubble solution is universal in the sense that 
the bubble action (and therefore the radius) indeed only depend on the surface tension, 
obtained from the one-dimensional domain wall problem, and the energy density splitting 
of the two vacua. The critical radius of the bubble, obtained by maximizing the action 
with respect to R, and the value of the action become 

Rc = 6 x t%' 5bubblc = w x \Ty~g ) ' (5 ' 4) 

Quantitatively, the applicability of the thin wall approximation is determined by the con- 
dition that the radius of the bubble is much larger than the width of the domain wall, a 
quantity of order one in the dimensionless units. Using the full theory value for the surface 
tension, this translates to <5 < 1, which is certainly satisfied for quarks with rrij > T, fij, as 
5 is then exponentially suppressed. 



6 Conclusions and outlook 

In the paper at hand, we have constructed a dimensionally reduced effective theory for two- 
color QCD, generalizing the earlier works of refs. [2, 11] by including in the consideration 
the effects of fundamental quarks of in principle arbitrary masses and chemical potentials 
on the dynamics of the Wilson line. In the process, we also identified and corrected an 
important algebraic error in the original papers, which led to qualitatively new understand- 
ing of the structure of center symmetric effective theories. In particular, we found that 
the tension and profile of a domain wall separating the two Z(2) minima cannot be used 
as matching conditions between the effective and full theories, but that instead become 
genuine predictions of ZQCD. 

The effective theory we constructed is formulated in terms of a 'coarse grained Wilson 
line' type variable Z defined by eq. (2.3), as well as the Lagrangian of eqs. (2.4)-(2.6). 
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The theory has by construction a notion of the center symmetry transformation of the full 
theory, and is invariant under it in the absence of dynamical quarks. The matching of 
the effective theory to the full one was performed by requiring that the former reproduce 
the correct long distance physics of the latter, a task most conveniently accomplished by 
demanding that the theory reduces to EQCD upon integration out of the momentum scale 
ttT. This we carried out explicitly in section 4, where the Sj parameters appearing in the 
effective theory Lagrangian were determined with the exception of one, for which a higher 
order computation is needed. 

Upon fixing its soft parameters, our effective theory becomes fully predictive, as the 
physics of the distance scales l/(gT) and larger is to a very good accuracy independent 
of the values of the hard parameters hi [11]. Thus, our theory is immediately amenable 
to nonperturbative lattice simulations, with which one may study two-color QCD over an 
extensive range of temperatures, quark chemical potentials and masses. 5 To this end, let us 
note three encouraging facts: First, as shown in refs. [27, 28], dimensional reduction does 
not as such require the temperature to be the largest energy scale in the system, but rather 
only that a modest scale hierarchy exist between the scale ttT and the electric screening 
mass, translating at high densities to a parametric condition g/j, < irT. Second, there is a 
priori no reason, why our effective theory could not be extended all the way to the limit of 
massless quarks, although it then has little relevance for the phase transition itself, which 
is driven by chiral dynamics. And third, it should be recalled that the only reason we have 
in this paper chosen to study two-color QCD and not the physical case of three colors is 
notational and computational simplicity. If the predictions of the effective theory turn out 
to match two-color lattice data well, then an obvious topic for further investigations will 
be the generalization of the present work to full three-color QCD. 

Finally, another interesting topic for future work would clearly be to consider the 
relation of our effective theory to the strong coupling effective actions derived in refs. [13, 
29-31]. It appears that these two approaches are strongly complementary in the sense that 
they approach the deconfinement transition from opposite directions; whether this can be 
used to gain more insight into the dynamics of the transition itself remains to be seen. 
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A EQCD parameters in the presence of massive fermions 

The onedoop QCD effective potential evaluated in a static background Aq field was first 
determined in refs. [14-16]. Rewriting it in a form invariant under global SU(2) symmetry 
and parameterizing the gauge field as Aq = a - a/2, the result becomes 

2 



^M^} Mil))' 



1 + 2e-0te*±/*) cos EM + e -2p(e jk ± H ) 
2T 



(A.l) 



j=l ± 

where (•) denotes the fractional part of a real number ((x) = x — [x\), /x,- the set of (flavor) 
quark number chemical potentials, and ejk = * k 2 + m 2 the dispersion relation of the jth 



quark flavor. Also, we used the shorthand notation \a\ = y/a • a. Expanding this expression 
in powers of a around zero and subtracting the contribution of the static modes (amounting 
to the term cubic in a), one may readily identify the EQCD parameters 

2 _2g 2 T 2 f d*k 



2 \ — " v f ^ ^ 71 f \ 

2g l^J (2^ /( ^ fe ' ^' 



(A.2) 

2g A T g 4 f d 3 k 
X = "3^ + ¥ ^ J J^f 1 '"^^ 

where f(x,(x) = [f(x + fi) + f(x — fJ>)]/2, the prime denotes differentiation with respect to 
x, and f{x) = l/(e^ x + 1) is the Fermi-Dirac distribution function. 

Compact as the above expressions (A.l) and (A.2) are, one can further evaluate the 
integrals over the quark momentum analytically in terms of the modified Bessel function 
of the second kind, K n . Expanding the logarithm in powers of fugacity and using some 
identities for the Bessel functions, one obtains the result 

i^i \ 2 / / i -*i \ \ 2 



4T 2 ^ ~ (_!)n |g| 

H 7T / m i / — x—K2{nprnj) coshlnpuj) cos — — - 

tt z J n z 2i 



(A.3) 



where the sum converges as long as //j < rrij for all quark flavors. Analogously, one derives 
by differentiation analytic expressions for the EQCD mass parameter and quartic coupling, 

2 g ^T 2 g 2 Nt °° 
m 2 x = — ^2 m 2 y^(-l)" K 2 (n/3mj) cosh(n(3fij), 

71 3=1 n=1 (A.4) 
~ X = ~t 2 ~ + ^^E^E^^^Mm.OcoshCn^). 

j=l n=l 
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When the quark mass is parametrically larger than both the temperature and the respective 
chemical potential, the infinite series in eq. (A. 4) can be replaced by its asymptotic form, 



On the other hand, for massless quarks at vanishing chemical potentials, the integrals in 
eq. (A. 2) are readily evaluated analytically and one finds m 2 = (2g 2 T 2 /3)[1 + (iVf/4)] and 
A = (2 5 4 T/3vr 2 )[l - (N { /8)], in agreement with ref. [6]. 

As the above infinite series containing Bessel and hyperbolic functions will appear 
frequently in our results, it is convenient to introduce a shorthand notation, 

JV f oo 

± - ^(/3my) 2 ^2(±l) n n e K 2 (n^m j ) cosh(n/%), (A.6) 



n=l 



in terms of which the EQCD parameters (A. 4) take the simple forms m 2 = (2g 2 T 2 /3)[1 — 
(3k^/27t 2 )] and A = (2<7 4 T/37r 2 )[l + (k^/8)]. 6 To introduce one final piece of notation, 
observe that in the presence of dynamical quarks, the potential of eq. (A. 3) has only one 
global minimum (up to periodicity) at \a\ = 0, while the point \a\ = 2irT/g corresponds to 
a mere local minimum. The most important quantity carrying information on the explicit 
Z(2) breaking due to dynamical quarks is thus the energy density difference between the 
two minima. It can be encoded in a single dimensionless parameter 

, _ Ves(9\S\ = 2nT) - V cS (g\d\ = 0) _ 4 + 
B Center symmetry for the SU(2) gauge group 

The zero component of the gauge field Aq (r) (dependence on the spatial coordinates does 
not play any role in what follows and is thus suppressed) transforms under the local trans- 
formation s(t) G SU(2) as Aq(t) = s(t)Aq(t)s(t)^ + ^s(r)d T s(Ty . Under the same 
transformation, the (untraced) Wilson line operator defined in eq. (3.2) transforms as 
U = s(j3)Qs(0y . Let us now write the gauge field as Aq(t) = a(r) • a/2, and the most gen- 
eral SU(2) gauge transformation as s(t) = exp[i9j(r)n(r) • a], where n(r) is a unit vector. 
In this representation, the gauge field transforms as 

a =n{n ■ a) + [a — n(n • a)] cos 2(p — {ft x a) sin 2ip+ 

H — \2nLp' + n' sin 2<p — 2{n x n') sin 2 ip\ , 
9 



Note that at any given time, the infinite sum can be replaced by the corresponding integral expression. 
This is in particular necessary for reasons of convergence, if rrij < fij for some quark flavor. 
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where the prime denotes a derivative with respect to r. Demanding that the gauge trans- 
formation preserves the periodicity of the gauge field, a((3) = 3(0), leads to the conditions 

ip{P) = p(0) + Ntt, <f/(J3) = (f/(0), n(/3) = n(0), n'(/3) = n'(0), (B.2) 

up to an overall minus sign, which only matters if we require the parameters ip, n to change 
continuously with r. In either case, the unitary matrix s(r) satisfies s(/3) = (—1) N s(0) for 
some integer N, which is precisely a transformation of the Z(2) center of the gauge group. 

Let us now specialize to the Polyakov gauge, in which Aq is diagonal and independent 
of r. Which gauge transformations from the local SU(2) group preserve this structure? Ob- 
viously, n must point in the third direction at all times. Preservation of time independence 
of a then results in the condition ipt = const, and the admissible gauge transformations 
take the form n = (0,0,1) and c/?(t) = tp(0) + Nttt/(3. As a consequence, the gauge 
field transforms by a mere overall shift, 03 — > 03 + INnT/g. It is worth emphasizing, 
though, that this conclusion only holds in the Polyakov gauge, as otherwise the vector 
function a(r) transforms in a rather complicated manner. In any case, since the nontrivial 
center transformations correspond to time-dependent s(r), while there is no time in the 
three-dimensional effective theory, this effective theory must be augmented with a suitable 
definition of the center symmetry. 



C One-loop effective potential of ZQCD 

The derivation of the one-loop effective potential of the theory defined by eq. (2.4) follows 
closely appendix A of ref. [11], and we will therefore merely write down the result here. 
In practice, we choose II to point in the third color direction, II a = |n|5 a 3, relying on 
the SU(2) invariance of the theory. The effective potential in a general renormalizable 
gauge then consists of the tree-level contribution, the gluon and ghost loops, the loop in 
the mixed XH3 sector, as well as a separate contribution from the IT 2 loops, 



Vcfi = — Vt rC c + VA+gh + Vsn 3 + Vrii, 2! 
#3 

y trcc = 61S 2 + 6 2 n 2 + ciS 4 + c 2 (n 2 ) 2 + c 3 x 2 ri 2 + ^s 3 + d 2 sn 2 , 



VA+gh 

Vm 2 



mi 



^ 3/2 ), 



^n 3 = -7^ 



6vr v " n A 12vr 

~(2& 2 + 4c 2 n 2 + 2c 3 s 2 + 2d 2 s + en 2 ) 3/2 , 

D7T 



(mlf 2 + (m 2 f' 2 



where 



m 2 ± =h + b 2 + 6ciS 2 + 6c 2 n 2 + c 3 (S 2 + n 2 ) + (3di + d 2 )Z± 



± 



(C.l) 



61 - 6 2 + 6 Cl S 2 - 6c 2 n 2 + c 3 (n 2 - S 2 ) + (3di - d 2 )sj + 4n 2 (2c 3 S + d 2 ) 2 . 

(C.2) 



Using the parametrization of the couplings (2.8) and assuming that the background fields 
S,II are of order 0(^3) or smaller, we immediately obtain an expansion for the effective 
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potential to next-to-leading order in the coupling, 



l 



5l L 2 



h 



1 /v2 



(£^ + n 2 ) + ^(s 2 + ri 2 ) 2 



+ £in2 + £2:(n2)2 + S3S 4 + + ^ EI1 



4 i S 4 V 3 s 5 v ff2 



2 4 
1 



HL (2 _ e 3/ 2) _^_ ( ^ +eri2 )3/ 2 __i_ ( ^, 



(C.3) 



where m\ = hi + 3/i 2 (E 2 + IT) and mi = h x + /i 2 (S 2 + IT) 
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